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We study solitary wave solutions of the higher order nonlinear Schrodinger equation for the 
propagation of short light pulses in an optical fiber. Using a scaling transformation we reduce 
the equation to a two-parameter canonical form. Solitary wave (1-soliton) solutions always exist 
provided easily met inequality constraints on the parameters in the equation are satisfied. Conditions 
for the existence of iV-soliton solutions (N > 2) are determined; when these conditions are met the 
equation becomes the modified KdV equation. A proper subset of these conditions meet the Painleve 
plausibility conditions for integrability. 
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The propagation of nonlinear waves in dispersive media is of great interest since nonlinear dispersive systems are 
ubiquitous in nature. Propagation of ultrashort light pulses in optical fibers is of particular interest because of the 
common expectation that solitary waves may be of extensive use in telecommunication and even revolutionarize it. 
The existence of solitary wave solutions implies perfect balance between nonlinearity and dispersion which usually 
requires rather specific conditions and cannot be established in general. The objective of the present paper is to study 
the conditions under which the existence of solitary waves is guaranteed for ultrashort pulses. 

The propagation of light pulses in fibers is well described by the higher order nonlinear Schrodinger equation 
(HONSE) a partial differential equation (PDE) whose right hand side includes the effects of group velocity 

dispersion, self-phase modulation, third order dispersion, self-steepening, and self-frequency shifting via stimulated 
Raman scattering, respectively: 



E z = i( ai Eu + a 2 \E\ 2 E) + a 3 E m + a 4 (\E\ 2 E) t + a 5 (\E\ 2 ) t E. 



(1) 



When the last three terms are omitted this propagation equation for the slowly varying envelope of the electric field, 
E, reduces to the nonlinear Schrodinger equations (NSE), which is integrable (meaning it not only admits TV-solitary 
wave solutions, but that the evolution of any initial condition is known in principle) ||-^|. We call these iV-solitary 
wave solutions iV-solitons, and mean by this that the solitary waves scatter elastically and asymptotically preserve 
their shape upon undergoing collisions, just like true solitons. However, for short duration pulses the last three terms 
are non-negligible and should be retained. In general, the presence or absence of solitary wave solutions depends on 
the coefficients a appearing in Eq. (|l|), and therefore, on the specific nonlinear and dispersive features of the medium. 
Here, we reduce the HONSE to a two-parameter equation and derive a general solitary wave (1-soliton) solution. We 
determine conditions when 7V-soliton solutions exist and display the solutions. We also study the Painleve plausibility 
conditions for integrability and show that these are only a proper subset of the conditions necessary for iV-soliton 
solutions to exist. 

We begin by scaling the HONSE, letting E = bi A, z = b 2 ( : and t = b 3 r. Substituting into the HONSE we obtain 

= i((& 2 ai/^)A TT + {b\b 2 )a 2 \A\ 2 A) + (b 2 a 3 /b 3 3 )A TTT + (b 2 b 2 a 4 /b 3 )(\A\ 2 A) T + (b 2 b 2 a 5 /b 3 )(\A\ 2 ) T A. (2) 

Choosing 61 = {a\/ {a^l)) 1 ^ 2 , b 2 = a 2 /a 3 , and 63 = a 3 /a\, we can set the coefficients of the first, second and third 
terms on the right hand side of Eq. (||) to unity, so that the HONSE becomes 

A c = i{A TT + \A\ 2 A) + A TTT + j!(\A\ 2 A) T + j 2 (\A\ 2 ) T A, (3) 

where 71 = b 2 b 2 a 4 /b 3 = a 4 a\/a 2 a 3 and 72 = b\b 2 a$/b 3 = a^ai/a 2 a 3 . 

A solution to Eq. (|^) of the form A(£, r) = y(r + (3Q exp[«(«;C — fir)], with y real, exists regardless of the values of 
the parameters 71 and 72 as can be easily demonstrated by substituting this form into Eq. (^5j) and equating the real 
and imaginary parts of the resulting equation. This procedure yields the following necessary and sufficient condition 
on Q, and an equation for k in terms of 17: 

n= 3 ?fr^ 3 ' ( ^ 1/3) < (4) 
d(7i + 72) 

k- n 3 + n 2 = (f3 + m 2 - 2ft) (1 - 3fi). (5) 
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The function y(£) (£ = r + /?C) satisfies the ordinary differential equation 



ya = (P + 3n 2 -2n)y-^±^y\ (6) 

whose solution is given generally in terms of doubly periodic elliptic functions. For zero energy (i.e., for 
Vl - (0 + 3ft 2 - 2tt)y 2 + (371 + 2-f 2 )y A /6 = 0) we find the solitary wave solution 

A(C, r) = ( ^^f^ ) ^h- 1 [(0 + 3ft 2 - 2ft) V 2 ( T + 0Q] 

x exp{i[((0 + 3ft 2 - 2ft)(l - 3ft) + ft 3 - ft 2 )C - Or]}, (7) 

provided (/3 + 3ft 2 — 2ft) > and 371 + 272 > 0. Thus, solitary wave solutions always exist (in contrast with what is 
implied in Ref. J?]]) provided 371 + 272 > 0. For ft = (i.e., for 371 + 272 = 3) the solitary wave solution reduces to 
A((, t) = 2 1/2 ?7Cosfr 1 (7/(T + rfQ) exp(i?7 2 C) where n = = rj 2 > 0. 

The case of ft = 1/3 is very special, not only because Eq. M) (which should be written in the form 371 + 272 = 
3(1 — 7ift)/(l — 3ft)) is not applicable, but also because Eq. (j^) is expected to be integrable for this case, as we shall 



show below using a Painleve analysis [|3-10|. The Painleve condition for integrability (see below) 

7i = 3,72 = -3/2, (8) 

(which differs from the result claimed in Ref. Q]) yields ft = 1/3 and k = —2/27, and the solitary wave solution takes 
the form 

A(C, r) = {0- 1/3) 1 / 2 cosh^p - l/?,) 1 ' 2 ^ + 0C)] ex P H(2C/27 + r/3)]. (9) 

It is of interest to compare the solitary wave solutions for different values of parameters 71 and 72. All solitary 
waves have intensity profiles of the form I((,t) — \A((,t)\ 2 = 7 s cosh- 2 [(T + 0Q/t s \. The solitary wave width r s 
and its intensity I s are related by the expression I s t 2 = 6/(371 + 272) (where 371 + 272 > for a solitary wave 
solution). For comparison, recall that I s t 2 — 2 for the NSE solitary wave (0:3 = = a$ = 0), so that for equal t s , 
jHONSE^jNSE _ 3/(3^ -|- 272). The width r s and (the negative of the inverse-velocity in the coordinate system 
moving with the group velocity of the light pulse, i.e., the solitary wave velocity v s = v g — 0~ x where v g is the group 
velocity) are related by the equation r s = (0 + 3ft 2 — 2ft) -1 / 2 where ft is given by Eq. 

As is well-known, existence of solitary wave solutions does not guarantee existence of 7V-soliton solutions with 
N > 1. In the absence of the inverse scattering solution for the HONSE, the Hirota method [O], based on the cutoff 
of the Pade-approximation [Q A = G/F, F being real, is often useful. Direct substitution of this representation 
into the PDE in Eq. (|) can be expedited using Hirota D operators |ll[] (defined by its operation on bilinear forms 
Dt{f ■ 9) '■= (d/dt — d I dt') f (t)g(t')\t = t>) and by separation of the linear part of the PDE to yield: 

(D c - iD 2 T - D 3 T )(G ■ F) =0, (10) 
(G ■ F)[-iD 2 T (F ■ F) + i(G* ■ G) + (71 + l2 )D T (G* ■ G)} 

+D T (G ■ F)[-3D 2 T (F ■ F) + (3 7i + 2 72 )(G* • G)] - 0. (11) 

The standard algorithm is to further substitute a 2-soliton solution in a power series in e of the form 

G = e(exp(#i) + cxp(0 2 )) + e 2 G (2) + . . . , (12) 
F = 1 + eF {1) + e 2 F (2) + . . . , (13) 
0i = p£ + qiT + <t> u *=1,2, (14) 

and require that the series truncate. 

The Hirota method is well adapted to the case when Eq. ([ll]) reduces to a bilinear equation or can be split into two 
or more independent (and consistent) bilinear equations. Unfortunately, further splitting of Eq. ( p"l| ) into two bilinear 
equations is impossible in the general case (the direct naive splitting of Eq. ( pd| ) into two equations corresponding 
to the two expressions in square brackets [p"3|]'ij] is incorrect, since it gives pi — p 2 and qi — q 2 for the 2-soliton 
substitution ( |l2| ) which is forbidden). Direct analysis of the multilinear equation ([ll]) is not successful either. We 
shall try further reduction by substituting 
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G = Gexp[z(/< - Or)] (15) 



into Eqs. ( |l0| ) and (|llj). Choosing n = 2 — 3 we retain the structure of Eq. ([L0[), while Eq (|ll|) takes the following 
form: 

(G ■ F)[-i(l - 30)£ 2 (F ■ F) + i(l - 7l )(G* ■ G) + ( 7l + -y 2 )D r (G* • G)] 

+D T (G ■ F)[-3D 2 T (F ■ F) + (3 7 i + 2 72 )(G* • G)] = 0. (16) 

Selecting O = 1/3 or O = I/71 eliminates the terms oc D 2 (F -F) or (G* -G) in the first part of Eq. (pit), respectively. 
Further simplification can be achieved in the particular case 71 = 3 when the choice = 1/3 leads to the following 
equation 

D T (G ■ F)[-3D 2 (F ■ F) + (9 + 2 72 )(G* • G)] + (G • F)(3 + 72)^r(G* • G) = 0. (17) 

Note that the first Hirota equation ([l0|) implies that any ./V-soliton solution which takes the Hirota form can be written 
as G/F = A("/2)G/F, where Q and F are independent of the parameter 72 (cf. p^pffl ) . Substituting this into ( |l7j ) 
shows that this is possible only if 72 = —3 or \n(Q* /Q) = const. The first choice corresponds to the case described in 
Rcf . [jll| . We shall consider here the second choice which leaves 72 arbitrary. 
Returning to Eqs. (||) we substitute 

A = £(£, t) exp[i( K C - Or)]. (18) 

which gives: 

£ c = -i(n + Q 2 - Q 3 )£ + (20 - 30 2 )£ T + i(l - 3fl)£ TT 

+i(l - njl)\£ \ 2 £ + £rrr + Jl(\£\ 2 £)r + l2£(\£\ 2 )r- (19) 

Choosing n = O 3 — O 2 and making the coordinate transformation T = r + (20 — 30 2 )£, eliminates the terms oc £ 
and £ T . Selecting = 1/3 or = l/ 7 i, eliminates terms oc £ TT or \£\ 2 £, respectively. In the particular case 7 i = 3, 
O = 1/3, one arrives at the following complex modified KdV equation: 

£ c = £ TTT + 3{\£\ 2 £) T + 12£(\£ | 2 )t, (20) 

which is more general than the equation considered in Ref. |l4| since it is for arbitrary 72. It is easy to see that, for 
the special case 71 = 3, Eq. (||) has an ./V-soliton solution of the form: 

A = B(C,T)exp[i(-2C/27-r/3 + ip)], T = t + C/3, V = const, (21) 

where B((,T) is the iV-soliton solution of the real modified KdV equation 

B c = Bttt + (9 + 2 l2 )B 2 B T , (22) 

and can be written as § 

B=i(^ W / /2 (Hr/f))T, (23) 

N N 

f= °p(Ew(<it+W2)+ /'•/'.- ( 24 ) 

li=0,l i=l l<i<j 

m = + qfC + vf\ M 3 = ln(*-* ) 2 . (25) 

qi + 1] 

This iV-soliton solution differs from that proposed in Ref. (which does not fulfill the corresponding nonlinear 
equation). Note that the existence of iV-soliton solutions does not in general imply integrability. However, the 
integrability of the complex modified KdV in the special case of 71 =3 and 72 = —3/2 has been demonstrated in 
Ref. @. 

Finally, we study the integrability of the HONSE, applying the Painleve analysis P pp[|. It is widely believed that 
possession of the Painleve feature is a sufficient criterion for integrability (see discussion in Refs. |l^,[l(|). The PDE 
in Eq. (3) can be analyzed to ascertain whether it is integrable by seeking a solution of the PDE in the Painleve 
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form A((,t) = (a(£,"r))~ CT J2m=o b m (0( a ((> r ))™\ where b m (() are analytic functions of ( in the neighborhood of a 
noncharacteristic movable singularity manifold defined by a(£, r) = t — /(C)- Following Ref. Jl6| we substitute into 
Eq. (3) the following Laurent series 

A = (T-T Q )-^J2 R m(C)(r-r a r, (26) 

771—0 

A* = (r - r )- ff £ S m(Q(r - r ) m , (27) 

m— 

(in the vicinity of a movable singular point r = to(C)- Inspection of the strongest singularity immediately gives a = 1 
and i?o<So = — 6/(371 + 272). The subsequent substitution of ( |26| ) and j2^ ) into Eq. (3) should allow identification of 
the other R m and S m leaving exactly 6 arbitrary functions (including tq(z)) undetermined. Analysis of resonances 
is carried out by making an auxiliary substitution of the form A = Rq(t — tq)^ 1 + R p (t — ro) p ~ 1 (and similarly for 
A*) into the PDE and retaining only linear terms in R p (cf. Ref. (|^1). This procedure shows that resonances occur 
ai p — —1,0, 3, 4, and 3 ± 2a/— 71/ (371 + 272). The resonances at p = — 1 and p = correspond to the arbitrariness 
of To and Sq/Rq, respectively. Requiring that the resonant indices p be integers one finds that either 71 = —272 
or 71 = —72. Further substitution of the complete power series of (^6|) and (|2?j) into Eq. (||) shows that there are a 
sufficient number of arbitrary functions if and only if 71 = 3 and 72 = —3/2. This corresponds to the case studied 
in Ref. |l4f |. Thus, the parameters that insure the Painleve feature are a proper subset of the parameters that insure 
existence of an 7V-soliton solution of the HONSE. 

To summarize, we have derived a 1-soliton solution of the general HONSE, without any constraints on its coefficients 
except for the weak inequality constraint 371+272 > 0. The 1-soliton width and intensity are related by the expression 
I s t^ = 6/(371+272), so the intensity decreases with increasing (371+272) for a given r s . We have shown that iV-soliton 
solutions exist when 71 = 3 and 72 is arbitrary, thus extending significantly the results of Ref. [jllLH • The Painleve 
plausibility condition for integrability, 71 = —272 = 3, has been shown to be a proper subset of the conditions for 
Af-soliton solutions, and consistent with the integrable case found in Ref. EJ, while integrability of our more general 
case is still unknown. Our analysis does not exclude existence of iV-soliton solutions nor does it disprove integrability 
for arbitrary 71 and 72. 
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